Global persistence exponent of the double-exchange model 
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We obtained the global persistence exponent 8 g for a continuous spin model on the simple cubic 
lattice with double-exchange interaction by using two different methods. First, we estimated the 
exponent 9 g by following the time evolution of probability P(t) that the order parameter of the 
model does not change its sign up to time t [P(t) ~ t~ 9g ]. Afterwards, that exponent was estimated 
through the scaling collapse of the universal function L 9gZ P(t) for different lattice sizes. Our results 
for both approaches are in very good agreement each other. 
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During the last few years, a great number of experi- 
mental and theoretical works have been published about 
manganites with perovskite structure, Ai- x B x Mn03, 
where A is a rare earth such as La, Nd, and Pr, B is a 
divalent element such as Sr, Ca, Ba, and Pb, and x is the 
concentration of B. Such attention is due to the redis- 
covery of colossal magnetoresistence (CMR) Ql H, H, Q|, 
an extremely large change in the resistivity when a mag- 
netic field is applied near of Curie temperature. Besides, 
these manganese oxides possesse metal-insulator (MI) 
transitions, as well as a rich variety of physical proper- 
ties and possible technological applications, for instance, 
magnetic sensors and memory technology. 

In these compounds, the MI transition associated with 
the ferromagnetic spin alignment has been widely ex- 
plained by the double-exchange (DE) mechanism 0, El 
0, Mi LJ- LLOl 7 m which the transfer of an itinerant e g elec- 
tron between the neighboring Mn ions through the O 2- 
ion results in a ferromagnetic interaction. Anderson e 
Hasegawa showed that the transference element is pro- 
portional to cos(0/2), where (f> is the angle between the 
ionic neighbouring spins. This result was recently con- 
firmed for layered manganites • Although this theory 
has also been succeeded in explaining qualitatively CMR, 
some authors have argued that it cannot alone provide a 
complete description of this phenomenon. They suggest 
that, in addition to the double-exchange, a complete un- 
derstanding of these materials should include strong elec- 
tron correlations [T^ . a stron g ele ctron-phonon interac- 
tion [l^, or coexisting phases |14j . One might therefore 
think that double-exchange alone cannot explain CMR 
in manganites |13| . but this remains an open question. 
What we know is that in the study of the manganites the 
double-exchange theory plays an important role, both in 
the study of CMR and in explaining the presence of a 
ferromagnetic state (for x ~ 0.3) in doped manganites, 
furnishing the basis for describing manganites with colos- 
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sal magnetoresistence. 

The Hamiltonian of a classical s pin model with double- 
exchange interactions is given by [15J 

n = -j ^i + Si -Sj, (i) 

<*,j> 

where indicates that the sum runs over all nearest- 
neighbor pairs of lattice site, J is the ferromagnetic cou- 
pling constant and the spin Si = (Sf, S? , S%) is a three- 
dimensional vector of unit length. 

The critical properties of the DE models have been 
intensively studied by using Monte Carlo (MC) simula- 
tions. In the equilibrium |15l llo. Il7l Il8|. the estimates 
for the static critical exponents indicate that this model 
belongs to the universality class of the classical Heisen- 
berg model [l^. Several experimental works about the 
critical properties of the doped perovskite manganites 
also support this assertion [gjj, El H2- Very recently 
the dynamic critical behavior of the DE model was stud- 
ied by using the short-time Monte Carlo simulations and 
estimates for the static critical exponents v and (3 and 
the dynamic critical exponents z and 9 were derived |23j| . 
That approach is based on the results of Janssen et al. 
|24| that showed that universality and scaling behavior 
are already present in systems since their early stages of 
the time evolution after quenching them from high tem- 
peratures to the critical one. By using renormalization 
group techniques they obtained for the fc-th moment of 
the magnetization, extended to systems of finite size |25l , 
the following scaling relation 

M (k) (t, e, L, mo) = b' kp/u M {k) {b- z t, b 1/u e, b^L, b x °m ), 

(2) 

where t is time, b is an arbitrary spatial rescaling fac- 
tor, e = (T — T c ) /T c is the reduced temperature and L 
is the linear size of the lattice. The exponents f3 and 
v are the equilibrium critical exponents associated with 
the order parameter and the correlation length, and z 
is the dynamical exponent (r ~ £ z where r is the time 
correlation). For a large lattice size L and small initial 
magnetization mo at the critical temperature (e = 0), the 
magnetization is governed by a new dynamic exponent 9, 
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M(t) ~ m t e (3) 

if we choose the scaling factor b = t l l z in Eq. J5J). 
This new critical index, independent of the previously 
known exponents, characterizes the so-called "critical ini- 
tial slip" , the anomalous increase of the magnetization 
when the system is quenched to the critical temperature 
T c . In the sequence, another new independent dynamic 
critical exponent was found by Majumdar et al. 26] to 
describe the behavior of the global persistence probability 
P(t) that the order parameter has not changed its sign 
up to time t. At criticality, P(t) is expected to decay 
algebraically as 

P(t)~ t - e °, (4) 

where 9 g is the global persistence exponent. Since then, 
the study of the persistence behavior have attracted an 
enormous interest, playing an important role in the study 
of systems far from equilibrium pfj Eli EM EM IM Isil 

m m h m m mm 

Since the time evolution of the order parameter is in 
general a non-Markovian process, the new critical expo- 
nent 6 g is independent of the usual exponents. However, 
as argued by Majumdar et al. j2(| , if the global order pa- 
rameter is described by a Markovian process, 9 g is not a 
new independent exponent, being related to other critical 
exponents by the equation 



trial orientation of the spin is randomly chosen and ac- 
cepted or rejected according to the probability e _/3 ^ B 
where E 1 (E) is the new (old) energy of the spin system, 
(3 = J/ksT c and ks is the Boltzmann constant. The 
estimates were obtained from five independent bins with 
5000 samples each for t up to 500 Monte Carlo sweeps. 

In the first method, we used the scaling relation given 
by Eq. J3J in order to obtain 9 g as a function of mo 
for several values of the initial magnetization. Here, it 
is necessary working with a precise and small value of 
the magnetization, mo ~ 0. Next, the estimate for 9 g is 
obtained extrapolating that series to the limit in which 
mo — > 0. 

In Fig. n]we show the decay of the global persistence 
probability for L — 60 and mo = 0.00125 in double-log 
scales. Fig. [5] exhibits the behavior of the exponent 9 g 
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The global persistence probability P(t) can be defined 

as 

t 

P(t) = 1-X>(*0. (6) 

t'=i 

where p(t') is the fraction of the samples that have 
changed their state for the first time at the instant t' . 

In this paper, we performed short-time Monte Carlo 
simulations to explore the scaling behavior of the global 
persistence probability P(t), for a classical ferromagnet 
with double-exchange interaction. The dynamical expo- 
nent 9 g that governs the behavior of P(t) at criticality is 
obtained by using two different approaches: the straight 
application of the power law behavior [see Eq. (JIJ] and 
by means of time series data collapse. To our knowledge, 
this is the first time that this exponent is calculated for a 
three dimensional model and continuous spin variables. 

In our simulations, we considered L x L x L (L = 20, 
25, 30, 35, 40, 50 and 60) simple cubic lattices with 
periodic boundary conditions. Simulations were done 
at critical temperature T c = 0.74515, in units of 
J/ks, where fee is the Boltzmann's constant. The update 
we used is local and follows the Metropolis algorithm, 
i.e., at each site of the lattice (during the simulation) a 
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FIG. 1: Time evolution of the global persistence probability 
Pit) for a lattice size L = 60 and mo == 0.00125. The error 
bars were calculated over 5 sets of 5000 samples. 

for m = 0.005, 0.0025 and 0.00125, as well as a linear 
fit that leads to the value 9 g = 0.336(9). 

In Table|Uwe present the estimates for 9 g in function of 
different initial magnetizations mo for other lattice sizes. 

Our results obtained through the linear extrapolations 
for m-o — > are presented in the last column. The es- 
timate of the exponent 9 g for odd lattice sizes (L = 25 
and L = 35) were obtained for different values of m . 
For L = 25 we used m = 0.0064, 0.0032 and 0.0016 
and the extrapolated value was 9 g = 0.337(5), whereas 
for L = 35 we used m = 240/(35 3 ), 120/(35 3 ), and 
60/(35 3 ) and the extrapolated value 9 g = 0.338(5). 

In the second method, we used the fact that the de- 
pendence of P(t) on the initial magnetization can be cast 
in the following finite-size scaling relation [2|J 

P(t) = t- e *f(t/L z ) = L- e °*f(t/L*), (7) 

where z is the dynamical exponent. Thus, the quan- 
tity L e « z P{t) is an universal function of the scaled time 
(t/L z ) and the wanted value of 9 g is that which fulfils 
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TABLE I: The global persistence exponent 9 g from the power law behavior for different initial magnetizations mo and even 
lattice sizes. 



L 


m = 0.005 


m = 0.0025 


m = 0.00125 


Extrapolated value 


20 


0.315(14) 


0.326(8) 


0.337(12) 


0.342(4) 


30 


0.315(5) 


0.322(12) 


0.335(9) 


0.339(5) 


40 


0.315(9) 


0.321(31) 


0.336(28) 


0.339(7) 


50 


0.319(10) 


0.325(11) 


0.336(14) 


0.338(6) 


60 


0.313(7) 


0.315(14) 


0.333(10) 


0.336(9) 



6 = 0.336(9) 
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FIG. 2: Dynamical exponent 6 g as a function of the initial 
magnetization mo for cubic lattices with L — 60. Each point 
represents an average over 5 sets of 5000 samples. 



FIG. 3: Scaling Collapse of the scaled persistence probability 
versus scaled time for Li = 50 and Li = 60. The curves were 
obtained from five independent bins of 5000 samples. 



that condition for different lattices. The best estimate 
for 9 g is found through the \ 2 test [39^. 

Unlike the first method used in this paper, in the col- 
lapse method the exponent 9 g is obtained without the 
need of careful preparation of the initial magnetization 
mo nor the limiting procedure. The only requirement is 
that (too) ~ where (•) is an average done over the sam- 
ples at t = 0. On the other hand, the collapse method 
demands the dynamical exponent z to be known before- 
hand. In this paper, we used the estimate obtained very 
recently for this exponent, z = 1.975(10) p3| . 

In Fig. |3| we show the collapse of the curves obtained 
for L — 50 and L = 60. The open circles show the 
collapse of the larger lattice rescaled in time. Our best 
estimate for 9 g for L2 = 60 and L\ = 50 is 

6 g = 0.335(9) (8) 

The estimates for other lattice sizes are shown in Table 
ITT1 These results are in very good agreement with the 
estimates obtained directly from the power law predicted 
in Eq. Q}. 

Using, for instance, the result of Eq. JSJ) and the es- 
timates of the exponents 9, z, (3, and v obtained in Ref. 
[23| . both for the largest lattice size (L — 60), we ver- 
ify through Eq. JSJ), the non-Markovian aspect of the 



TABLE II: The global persistence exponent g for the best 
collapse for the double-exchange model. 



L 2 1— » 


U 


o 3 


60 1 — * 


50 


0.335(9) 


60 1 — ► 


40 


0.333(13) 


60 1 — y 


30 


0.329(14) 


50 1 — > 


35 


0.337(11) 


40 1 — ► 


25 


0.330(12) 


40 1 — y 


20 


0.332(12) 



phenomenon we are deal with, since 
6 g = 0.335(9) and a g = -B + ~ — = 0.026(17). (9) 

Thus, the global persistence exponent is also in this 
case independent of other critical exponents but the dif- 
ference between our estimate for 9 g and the value ob- 
tained from Eq. J^J is greater than that observed when 
discrete spin models were analyzed (see Table ITTT1) . 

In summary, we have performed short-time Monte 
Carlo simulations in order to investigate the scaling be- 
havior of the persistence probability P(t) for a three di- 
mensional system with double-exchange interaction. The 
dynamic critical exponent 9 g that governs the behavior 
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TABLE III: The exponents 9 g and a g for several models. 



Models 


6 9 


a a 


Ising model [28] 


0.236(3) 


0.212(2) 


Three-state Potts model [28] 


0.350(8) 


0.324(3) 


Blume-Capel model [30] 


1.080(4) 


0.904(21) 


DE model (see Eq.@) 


0.335(9) 


0.026(17) 



of Pit) at criticality was estimated by using two different 
approaches: the straight application of the power law be- 
havior P(t) ~ t~ 0g and the collapse method for the uni- 
versal function L 9gZ P{t). The results are consistent with 
the expected non-Markovian character of the process. 

This work was supported by the Brazilian agencies 
CAPES and CNPq. 
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